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Abstract 

In this paper, we describe the line Dirac delta function of a curve in three-dimensional space 
in terms of the distance function to the curve. Its extension to level set formulation and plane 
curves are also developed. 
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1 Introduction 

A line source or sink is often encountered in electromagnetic and biological problems when a long 
and thin structure is involved. For instance, a distribution of electric charge on a curve can be 
represented by a line source pQ. A flexible structure such as an elastic fiber or filament is modeled 
as a line source in Immersed Boundary method to address the mechanical interactions between the 
structure and biofluid A very thin blood capillary in a bulk tumor can be treated as a line 
source when modeling a growth factor's reaction and diffusion [8]. In these models, the line source 
is considered as a simple smooth curve S, embedded in the three-dimensional (3-D) Euclidean space 
IR 3 or its generalization. 

Denote the arc length parameter by s and the corresponding spatial point on this curve by x(s). 
The line delta function 5^, associated with the curve £ is defined as a distribution such that for any 
test function /(x) G C^°(M 3 ) (i.e., infinitely differentiable and compactly supported), 

/ fc(x)./(x)dx= / f(x(s))ds. (1) 

In Cartesian coordinates, by formally switching the order of integrations (i.e. applying "Fubini's 
Theorem"), it can be understood as 5s( x ) = f% 5 3D (x — x(s)) ds, Vx 6 IR 3 , where <5 3i:) (x) = 
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5(x)5(y)5(z) for x = (x,y,z), where S(x) is the one-dimensional (1-D) Dirac delta function satisfy- 
ing 6(t)g(t)dt = g(0) for any g G C C °°(M). 

The purpose of this work is to provide an alternative representation of the line delta function 
Ss for a fairly general class of spatial curves sufficient for most practical applications, which is 

where p(x) is the distance function to the curve. This 5 function is under the convention that 
/ °° 5(t)g(t)dt = 5 (0) for any g G C C °°(R). 

The main challenge is to how to make sense of and calculate J K3 ^^/(x)dx for any / G C^°(IR 3 ). 
For this, one has to choose the most reasonable and practical coordinate system. Indeed, the appar- 
ently singular term - on the right hand side comes naturally from the local cylindrical coordinates 
in a tubular neighbourhood near the curve. If the curve £ is simple and smooth, then a tubular 
neighbourhood can be readily chosen where the distance function works as the radius and will 
be cancelled out in coordinate transformation. This will be done in Section 2. In general, if the 
curve has non-smooth points, we have to make the meaning of more clear, which will be done 
in Section 3. As application, the formula ([2]) works for very complicated structure composed by 
piecewise smooth curves such as the newly formed blood capillary network in growing tumor, as 
shown in Fig.[TJ 




Figure 1: The highly irregular and tortuous blood vessel capillaries in a xenotransplanted U87 
human glioblastoma multiforme tumor (upper part) in a mouse brain. The size of tissue shown in 
this figure is 2.6mm by 2mm. This picture is taken from [6] with permission. 



Another important application of this work is in the level set method. The delta function of a 
curve or a surface represented by a level set function is crucial in the level set method [1], but the 
form for a spatial curve is not previously available. The formula (|2|) provides the first representation 
with the distance function, whose generalization to an arbitrary level set function will be given in 
Section 4. The extension to plane curves and general submanifolds with co-dimension 1 and 2 are 
also discussed in Section 4. 

2 Case I: simple smooth curves embedded in R 3 

The essential tool in our calculation of |^ is the knowledge of a neighbourhood around the curve. 
The standard tubular neighbourhood theorem is for compact smooth submanifolds of M. k (c.f. 
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[2]). Therefore, when it is applied to a curve, the curve needs to be a submanifold of R 3 and is 
diffeomorphic to a circle. However, because of the local nature of our problem, it can be readily 
applied to any simple smooth curve which is a submanifold of R 3 . More detail is as follows. 

Let X be a simple smooth curve with the arc length parameter s, any point on the curve denoted 
as x(s). Look at a finite piece of it with s G [sj, S2], for simplicity, still denoted as X. It is a closed 
set in R 3 with the two endpoints x(si) and x(s2). Denote X as the interior part by removing 
the two endpoints of the curve X. On each point x(s) G X, denote the tangent space as T X ( S )X 
which is exerted by Jj. Define the normal plane as N x f s \'S = {z G R 3 |z _L T x( - Sl )E}. On the two 
endpoints x(si) and x(s2), define tangent and normal planes as the limits of those of points in X. 
Let B(x(si),e) be the open ball centered at x(si) with radius e, and B out (x(si), e) be the open 
semi-ball of -B(x(si), e) cut by the normal plane iV^^X which is disjoint with X. Similarly we can 
define B out (x(s 2 ) , e) . 

Define the tubular neighbourhood 

Tub e (X) = {x G R 3 \dist(x,Z) < e}\{B out (x( Sl ),e)UB out (x(s 2 ),e)}, 

which is the set of all points within distance e to X but excluding the two outside semi-balls at the 
two ends of the curve. Note A r x ( Sl )S and iV x (,, 2 )X are two bounding normal planes of this tube (see 
Fig. [2}. Also define 

5(2, e) = {(x,z) G X x R 3 |z G iV x{s) X, ||z|| < e}. 

Let 8 : H(X, e) — > R 3 be given by #(x, z) = x + z. The following is a standard result in differential 
geometry. 

Lemma 1. There exists e > such that the map 9 is a diffeomorphism from H(X,e) onto Tub e Y>. 

Remark 1. A prominent feature of this tubular neighbourhood is that\/y G Tub 6 Y>, dist(y, X) = \ \z\\ 
where the vector z is from the unique decomposition y = x + z, (x, z) G H(X,e). Also notice that 
for our curve of infinite length, one might not have a uniform e for the whole curve, but this is not 
going to give us any trouble in this work because only pieces of finite length will be considered in 
the proof of Theorem [7J 

Define the distance function to X as p(x) = inf yg £ (i(x, y), where <i(x, y) is the standard Eu- 
clidean distance. /?(x) = implies x G X because the curve X is embedded in R 3 (i.e. being a 
submanifold). 

Theorem 1. Let X be a simple smooth curve embedded in R 3 , then 

Proof. Around any point on X, we have a tubular neighbourhood Tub e X satisfying the condition 
in Lemma [H First, a coordinate system will be constructed in Tub e X as follows. On the curve X, 
^ is a unit vector field. Because the curve is smooth, there exist two other smooth unit vector 
fields U and V such that U, V, and form an orthonormal system and satisfy the right hand rule. 
On each normal plane -/V x (s)X, U and V provide a R 2 coordinate system, (u,v). Denote any point 
on this plane as x(it, v, s) (see Fig. [2]). According to RemarkHJ p(x(ii, v, s)) = \Jv? + v 2 . 

Because {U,V, ^} is a smooth frame along the curve, the map from (u, v, s) to (x, y, z) is 
smooth. The Jacobian between (x, y, z) and (u, v, s) is the determinant of the transition matrix 
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Figure 2: A tubular neighbourhood Tub e E of a simple smooth curve E in 3-D space bounded by 
two normal planes iV x ( sl )£ and iV x ( S2 )E. 



between the M 3 standard basis J^, ^} and J^, J^} which is just {[/, V, J^} when restricted 
on the curve. Although is not explicit away from the curve, the transition matrix has to be 
in SO(3) on the curve because it is between two orthonormal bases with the same (right-handed) 
orientation. So the determinant is 1 on the curve. As the Jacobian is smooth in the tubular 
neighbourhood Tub e £, it's 1 + 0(Vu 2 + v 2 ) in sight of the Taylor series. Thus the Jacobian must 
be nonzero in a sufficiently small neighbourhood, say, Tub ei £,0 < ei < e. For simplicity, it is still 
denoted as Tub e £. Therefore, (u,v,s) is a coordinate system in Tub e E. 

When converted to the polar version (p, 9, s), since det ( a^'g'g j ) = Pi ^ i s eas y to see 

For any test function /(x) G C^°(M 3 ), J* R3 |^/(x)dx only depends on the situation near E 
as 5(p) is supported on £. Since / is compactly supported, we can assume that it is supported 
in a ball B. Because B n £ is compact, there exist finitely many local tubular neighbourhoods 
centered at involved pieces of £ covering B n E, denoted as {Tj}i 6 /. We can also make sure that 
{Ti} with Tj <s Tj also covers BnS. Now a stardard partition of unity construction provides us 
with hi G C£°(Ti) such that X^ie/ ^« = 1 when restricted to Ujg/Tj. Finally, we can carry out the 
following computation. 



2vrp 



E / L I ^f(P,0,s)h l (p,e,s)(p + O(p 2 ))dpd9ds 



s J0 J p 



E / LI 6 -^-f(P,0,s)h l (p,9,s)(l + O(p))dpd6ds 



i€l 



Is J0 J p 
<-2- 



iei 

In sight of the definition in Eq. (p}, the theorem is proven. 



E / / ^f(0,9,s)hi(0,e, S )d9ds 

/(x( S )) E ^i(x(*))da = ^ /(x( a ))cfe = | /(x(*))da. 
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□ 



Remark 2. The choice of coordinates does not have to be fixed and is actually flexible. However, 
we do need to start with a Euclidean orthonormal basis along £ with one vector in the direction of 
the curve. In short, the delta function would function on the normal plane of the curve with the 
Euclidean structure induced from R 3 . 

3 Case II: general curves embedded in K 3 

With a proper understanding of the delta function, the previous argument can be adjusted to a 
general class of curves which are topological graphs with smooth edges. The graph can have infinite 
vertices with finite edges between any two of them and the edges need not be straight (i.e. being 
topological). We only need to require that there is no local "accumulation" happening, which is in 
general what people think of graphs. 

Let £ be such a curve. The whole curve £ may be closed or open, may have non-smooth points 
and self- intersections, and may have more than one connected components (See Fig. [3]). Denote the 




Figure 3: This curve has two connected components: the left component is non-closed, has non- 
smooth points at A and B, and crosses itself at C; the right component is closed and smooth. 



collection of all non-smooth points as Py, = {Pj,j £ T}, where T is a finite or countable index set. 
To handle these non-smooth points, we remove from £ a small neighbourhood in M 3 around each 
such point such that the boundary of the neighbourhood is perpendicular to the curve at the point 
of intersection. Denote the union of these neighbourhoods by U e and shrink U e to Py as e — > 0. 
Let Xe be the characteristic function of M 3 \ U e , i.e., \ e = over U e and 1 otherwise. 
The following definition is what we want: 

S( P ) _ ^ Xe • 5{p) (4) 



Notice that Xe 27Tp P can be understood in the same way as before because a neighbourhood of the 



2-7T/9 e->o 2irp 

2np 

non-smooth points has essentially been removed. The limit is in the weak sense, i.e. in the sense 
of current, with the existence justified by the following computation. For any / G C^°(M 3 ), 

lim / *^M/( x )dx = hm / Xt ■ f(x(s))ds = [ /(x(a))d 8 

The first equality makes use of the same kind of computation as that in the previous section, i.e., 
using a partition of unity construction to reduce to local tubular neighbourhoods with cylindrical 
coordinates. This is done for each e and the local tubular neighbourhoods only need to cover T,\U e . 
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This also shows the independence of the choice (on Xe) i n the construction. Simply speaking, 
there is no contribution from Ps i n integration. In summary, we have proven the following main 
theorem. 



Theorem 2. Let E be a topological graph in M with smooth edges and p be the distance function 

S(p) 

2wp ' 



to E. Then the definition @ makes sense and o~e(x) ~ 



4 Extensions: level set function, plane curves, and more 

Let be a non-negative level set function of E, that is, E = {x | 0(x) = 0}. Assume (ft is a function, 
4> = 4>{p,9,s), locally near E in the coordinate system used before, and ^ ^ around E. Since 

4> is non-negative, we have actually assumed §^ > 0. This implies (4>, 9, s) can be viewed as a 
generalized cylindrical coordinate system locally near E, and the Jacobian between it and (p, 9, s) 
is g(p'g'g| = Therefore, we have the following result. 

Theorem 3. Let E be a topological graph in M 3 with smooth edges, p be the distance function to 
E, and (/>(x) be a non-negative level set function o/E as above, then <fe(x) = 2t ^ p . 

The key is that dpdOds = d( ^/Q p d(f)d9ds. For non-smooth points, one can have similar discussion 
as in Section 3. 

Analogous to Theorem 2 and 3, we can obtain the following results for plane curves. 

Theorem 4. Let E be a topological graph in M. 2 with smooth edges, p be the distance function to 
E, and </>(x) be a non-negative level set function o/E, with ^ is uniformly bounded away from 0, 
then «J s (x) = 5(p) = 5((f))\V4>\. 

Proof. We only need to prove the level set formulation for a simple smooth curve E embedded 
in M?. We use the 2-D coordinates (p, s) in the tubular neighourhood U around E. Over E, 
= PWp *75s = P^p as §s = over E, where p and s are orthonormal vectors in the direction 
of §- p and ^, respectively. Then over E (i.e. {x | p(x) = 0} = {x | 0(x) = 0}), |V0| = f£ > 0. 
Take any test function / 6 C C °°(R 2 ), 

f(x)6(<j>)\V<f>\dx = I /(x)*(0)|V0|Jdpds 
'u 

/(x)J(0)|V0|J^^#d S 
f(x(s))ds, 

s 

where J is the Jacobian between the Euclidean coordinate and (p,s), which is equal to 1 along E 
because and ^ form an orthonormal basis there. 

□ 



Remark 3. For Theorem 3, is not well-defined as a function over E in general because of the 
presence of 9, and should be understood in the integration sense. Because we do not have 9 for a 
neighbourhood of a curve in R 2 , this is not an issue for Theorem 4- 
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Remark 4. In the 3-D case, a succinct formula like -^j|V0| does not hold in general. Indeed, in 
3-D cylindrical coordinates, |V</>| ~ yj (j>p + ^ near S. Therefore, when approaching S, ^ 

cannot be replaced by |V</>| unless (fte/p — > 0, which puts restriction on this level set function <f>. 

Remark 5. The study here can be extended to the general case of a submanifold, S, of co- dimension 
2 or 1 in W 1 with n > 3 (or any smooth Riemannian manifold in general). That is, 5^, = for 
a submanifold E with co-dimension 2, while 5y, = S(p) = <5(0)|V0| i/S is o/ co-dimension 1. 

Indeed, the form 5s = ^(0)|V0| for smooth closed curves in M 2 and surfaces in M 3 was first 
derived in [3] and now is widely used in the level set method. Hence, the results presented in this 
work can be regarded as extensions of [3]. 
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